arXiv: 1504.04783vl [quant-ph] 19 Apr 2015 


Prospects for observing dynamical and anti- dynamical Casimir effects in circuit QED 

due to fast modulation of qubit parameters 

D S Veloso^ and A V Dodonov^’^ 

^ Institute of Physics, University of Brasilia, 70910-900, Brasilia, Federal District, Brazil and 
^ International Center for Condensed Matter Physics, 

University of Brasilia, 70910-900, Brasilia, Federal District, Brazil 

We consider the nonstationary circuit QED architecture, where a single artificial two-level atom 
interacts with a cavity field mode under external modulation of one or more system parameters. 

Two different approaches are employed to study the effects of Markovian dissipation on modulation- 
induced transitions between the atom-field dressed states: the standard master equation of Quantum 
Optics and the recently formulated dressed-picture master equation. We estimate the associated 
transition rates and show that photon generation from vacuum (“dynamical Casimir effect”, DCE) 
and coherent photon annihilation from nonvacuum states (“Anti-DCE”) are possible with the current 
state-of-the-art parameters. 

PACS numbers: 42.50.Pq, 42.50.Ct, 42.50.Hz, 32.80-t, 03.65.Yz 


I. INTRODUCTION 


The subject of photon generation in nonstationary 
circuit QED, where mesoscopic “artificial atoms” (con¬ 
structed with Josephson junctions [II, i) interact with 
the Electromagnetic field confined in superconducting 
stripline resonators under nonstationary conditions Q, 
has been studied theoretically for nearly ten years. Tra¬ 
ditionally one considered the externalljyjrescribed mod¬ 
ulation of either the cavi ty freq uency |4l-[l^. the atom- 
field coupling strength [la - lr^ or the atomic transi¬ 
tion frequency [l8l - l23l | , while recently more sophisticated 
schemes involving multi-level atoms, several qubits or 
coupled cavities were examined [ifl, [23 - 1^ . Such a 
unique possibility of in situ control of the fundamental 
parameters of the Hamiltonian became possible thanks to 
the highly controllable solid state environment in which 
the system properties can be manipulated by external 
electric and magnetic fields [29l - [^ . Recently it was 
shown that modulation of any system parameter gives 
rise to similar effects of roughly the same order of mag¬ 
nitude, and simultaneous modulation of different param¬ 
eters with the same frequency can increase or decrease 
the associated transition rate depending on the relative 
phases of modulations [ll, Moreover, “multi-tone 

modulations” comprising two or more harmonic functions 
of time can lead to new types of effective interactions 


However, actual circuit QED architectures suffer from 
unavoidable dissipation effects, in particular, the damp¬ 
ing of the cavity field, the atomic relaxation and the pure 
atomic dephasing [36|. Although some works on photon 
generation in nonstationaiw circuit QED analyzed the 
effects of losses i ill Hi, ill, 113 , the majority of stud¬ 
ies investigated only the unitary dynamics for different 
regimes of parameters and modulation shapes. Now that 
the general analytical description in absence of losses has 
been formulated in a closed form (under a series of ap¬ 
proximations) [ 33 , [13 , the question of utmost practical 


interest is how dissipation affects the modulation-induced 
phenomena and whether they can be implemented with 
present or near-future technology. 

So the goal of this work is to unveil which nonsta¬ 
tionary phenomena can be implemented experimentally 
with the current technology [38l - l40l | and comprehend how 
different dissipation channels affect the time evolution. 
We concentrate on the single-qubit, single-mode circuit 
QED setup and take into account the three aforemen¬ 
tioned dissipation mechanisms using the Markovian and 
zero-temperature approximations. The master equation 
appropriate to this case, which takes into account the 
qubit-resonator coupling, was deduced microscopically in 
[Sq . In this so called “dressed-picture master equation” 
(DPME) the dissipation rates depend on the spectrum of 
noise evaluated at the frequency matching the transition 
between the atom-field eigenstates (dressed states), so a 
prior knowledge of the bath structure is required. Be¬ 
sides, the analytical solution is complicated in majority 
of cases, so we carry out exact numerical simulations by 
assuming that the noise has a flat spectrum at positive 
frequencies and is zero for the negative ones. The details 
of the underlying mathematical formalism can be found 
in section im 

We consider several regimes of photon generation and 
annihilation due to modulation of the atomic transition 
frequency and the atom-field coupling strength. Eirst we 
study the generation of two excitations from vacuum in 
the resonant regime (section IHI Al) and the creation of 
one photon and one atomic excitation in the dispersive 
regime {Anti-Jaynes-Cummings behavior, section [111 Bl) 
for single-tone modulations. Then we analyze the re¬ 
cently discovered Anti-dynamical Casimir effect (Anti- 
DCE, section PlI Cl) , whereby two excitations are coher¬ 
ently annihilated due to the modulation of qubit param¬ 
eters in the dispersive regime [13 • We show that the 
above phenomena can be enhanced by using two-tone 
modulations fsection llVI) . for which a total of four ex¬ 
citations can be generated from vacuum or annihilated 
from a known initial state. Quantum states with even 








2 


more excitations can be generated from vacuum in the 
dispersive regime via the dynamical Casimir effect (DCE, 
section 0, when the system parameters are modulated 
with frequency roughly equal to twice the cavity unper¬ 
turbed frequency. 

It is shown that all the above phenomena can be im¬ 
plemented with current technology in dissipative circuit 
QED, although DCE and Anti-DCE require state-of-the- 
art dissipation rates and precise tuning of the modulation 
frequency. Moreover, we demonstrate numerically that 
for time intervals of interest the predictions of DPME 
are almost indistinguishable from the ones obtained for 
the “standard master equation” (SME) of Quantum Op¬ 
tics [ 411 , 113 , which is much simpler to handle and whose 
approximate solution is deduced here in the long-time 
limit. 


II. MATHEMATICAL FORMALISM 

The atom-held interaction is described by the Rabi 
Hamiltonian [2l|, |4^ and for generality we take into 

account the parametric amplihcation term due to even¬ 
tual time-modulation of the cavity parameters [d^. The 
total Hamiltonian reads (we set = 1) 

H = u;n+^\e){e\+g{d+d'’'){<7++a-)+ix{d''''^ — a^), (1) 

where d and d^ are cavity annihilation and creation op¬ 
erators and n = d^d is the photon number operator; 
( 7 + = \e){g\ and d- = \g){e\ are the atomic ladder opera¬ 
tors, where \g) {|e)) denotes the atomic ground (excited) 
state, uj is the cavity frequency, D is the atomic transition 
frequency, g is the atom-held coupling strength and x is 
the squeezing coefficient related to the parametric am¬ 
plihcation process. The time-independent part of x rnay 
appear due to the terms proportional to the square of 
the vector potential [d^, while the time-dependent part 
is due to the time-modulation of the cavity frequency and 
in the simplest case of DCE reads x = (da;)“^da;/dt. 

In this work we follow the convention of papers [s^, 
and suppose that all the system parameters can be per¬ 
turbed simultaneously by external modulations of the 
form X = Xq + ex sin( 77 ('^^t -I- Here 

X = {a;,D,g,x}, Xq and ex > 0 are the corresponding 
bare values and modulation depths, and the sum runs 
over all the fast modulation frequencies > wq. We 
consider small modulation depths given by inequalities 

5 where is the max- 

lb 

imum number of excitations. The parameters ^ 0 
(7) 

and (px relative weights and phase constants 

corresponding to the modulation of parameter X at fre¬ 
quency . To shorten the notation we introduce the 
complex modulation depth e^x — ^xw^x ®xp(*t!>x^) that 
incorporates both the weight and the phase for the mod¬ 
ulation frequency In our numerical examples the 

phases cfx always be set to zero and Wx = 1 for 


single-tone modulations with frequency . Whenever 
more than one parameter are modulated simultaneously 
we shall call the process “multi-modulation”. 

We assume weak atom-field coupling, gQ^/nffffff <C wq, 
and xo = 0 , so one can describe the dynamics in the basis 
of eigenstates of the bare Jaynes-Cummings Hamiltonian 
(JCH) Hjc = woh -I- Do|e)(e| -I- 5 o(d(J+ -I- [HI, 
These states, known as dressed states, read 

|<bn,s) = s„,s| 5 ,n)-|-c„, 5 |e,ri-l) , ri > 0 ,5 = ± , ( 2 ) 

where n is the total number of system excitations and 
the index S labels different eigenstates with same n. 
We define formally |</?o,-i-) = 0 introduce the nota- 
tion Sm,— — Om,-t- — eosO^ryi, SYn,-\- — — sin^TT^, 

where Oq = 0, 0„>o = arctan[(A_ -|- fin)/{2go^/n)] and 

Pn = The eigenenergies are Aq,- =0 and 

Xn,s = Won —A_/2-|-5/3ji/2, where A_ = wq —Dg is the 
bare detuning (we assume |A_| ^ uiq). 

The approximate unitar y d ynamics for the Hamilto¬ 
nian o was solyed in [^ for arbitrary small-depth 
modulations. In the interaction picture defined by the 
unitary transformation Ut = exjp(—itHjc) the effectiye 
Hamiltonian reads 

H = - 

m,S,T j 

'x\‘Pm,T){‘Pm+2,s\+h.C. , (3) 

where m > 0, Sj'T = ± and the time-independent 
coefficients given throughout the paper 

(the interaction-picture density operator is p = TjlpUt). 
By choosing 77M = Xm+2,s — Xm,t one resonantly cou¬ 
ples the dressed states |<pm,t) \y^M+ 2 ,s)- If in ad¬ 
dition |Am-i- 2 ,s - Am^r - ??m| > \^m]r.,s\ val¬ 

ues of {m,S,T}, all other terms can be neglected un¬ 
der the Rotating Wave Approximation (RWA) (sj, 
and one obtains a simplified time-independent Hamilto¬ 
nian. Actually, the eigenenergies Xn,s rnust be corrected 
by the so called “intrinsic frequency shifts” of the or¬ 
der O(po/wo, Ex/wo) with X = {uj,Q,g,x}- These cor¬ 
rections do not alter significantly the mathematical for¬ 
malism, so we neglect them in the analytical derivations 
below; however, they will be incorporated into our nu¬ 
merical simulations to achieve exact resonances. 

In the presence of dissipation the dynamics must be 
described by the master equation for the density operator 

dp/dt =—i[H, p] +/fp, (4) 

where >6 is the Liouvillian superoperator. We shall use 
two different Markovian kernels to estimate the effects 
of dissipation. The quantum optical “standard master 
equation” (SME) at zero temperature is [4lj | 

^SME* = KD[d] • +xV[a-\ • +^V[aP\* , ( 5 ) 

where 'D[0]p = \{20p0^ — O^Op — pO'^O) is the Lind- 
bladian superoperator and az = |e)(e| — \g){g\- The con¬ 
stant parameters k, 7 and denote the cavity damping. 
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qubit relaxation and the qubit pure dephasing rates, re¬ 
spectively. 

The second approach is the “dressed-picture master 
equation” (DPME) at zero temperature developed in [s^ 


Ldr* = 




l,k^l 


+ Y, (rjf+ r!f)p[|z)(fc|]., 

l,k>l 


( 6 ) 


where we use the shorthand notation |^) to denote the 
JCH dressed states |<Pn,s), I increasing with the en¬ 
ergy Xn,s- The parameters of equation (jb)) are de¬ 
fined as’ch' = [70(O)/2]i/2a“, = 70(Afc,)l4"l V 2 , 

and ■ Here K^iy), 

"fiiy) and 70 (!z) are the dissipation rates correspond¬ 
ing to the resonator and qubit dampings and dephasing 
noise spectral densities at frequency v; we also defined 
Aki = Xk - Xi, = {l\az\k), = (l|(a -I- a^)\k) and 

= (^|(i5'+ + ( 7 _)|fc). We do not consider a specific 
model for the reservoirs and make the simplest assump¬ 
tion that for V > 0 the dissipation rates are constant and 
equal to the corresponding rates of SME, while for < 0 
they are zero. 

We shall solve both master equations exactly via nu¬ 
merical integration of the differential equations for the 
density matrix elements. We assume the usual value 
ijJqI'I'k = 8 GHz for the cavity frequency, high but fea¬ 
sible value go = 5 x 10 “^wq for the qubit-field coupling 
strength (within the weak coupling regime), moderate 
range of detuning |A_| < 8go and the state-of-the-art 
dissipation rates 10“®go [38l - l40l| . We 

presume the tentative value |T| ^ 10“^ for the collective 
relative modulation depth [T ^ see equa¬ 

tions (froi) . (HH), (l?71) . (itSl) below], small enough to fulfill 
the range of validity of the Hamiltonian ([3]) yet sufficient 
to implement the phenomena of interest. 


A. Simplifications for the standard master equation 

The interaction-picture density operator obeys the 
master equation 

^ = -ilH, p] + KV[d]p + 'yT>[d-]p + ^V[d^]p, (7) 

where O = UjOUf The natural basis to expand p con¬ 
sists of the JCH dressed states 

OO 

p(f) = XI X! PnJn{i)Wn,s){^m,T\ ■ (8) 

n,m—0 ih 

By substituting equation (|5]) into © one can obtain 
the differential equations for the density matrix elements 
Pn'mit)- For example, for the pure atomic dephasing we 


get [using p„;+ = (p +’_„)1 




N,M 


dt 


70 

2 1 


(cos 20m cos 20n - 1 ) Pn% 
-ttPM 20M cos 20Nptj'Ji 

sin 20N cos 20m iPM~N )* 
^g*t(/3A.-/3M) 20M sin 20Np--j^ 


—e 


(9) 


dpN 


N,M 


dt 


70 

2 L 


(cos26>MCOs26'Ar - (10) 


+e 


-itpN 
JtpM , 


sin 20 N COS 20 Mptf\M 


sin 20 m cos20NipJi~j^) 


+e 


-U^Pn-Pm) 


-t,— '1* 
m , n ) 

sin20 m sm 20 Np ^ r'ti 


dPN’,M 

dt 


70 

2 1 


{cos20mcos20n+ 1)p%’m ( H ) 


+e 


itPu 


sin 20 m cos20mP%'% 

sin 20 N COS 20 mp]v’m 


- e < P ^+ P ^) sin20 m sin20x(/5]F]v)‘ 


When one considers the pair of equations 
j^A{t) = -qe^*^B{t) , j^B{t) = -qe-^*^A{t) (12) 

with constant parameters q and re, the solution is 

^ ^ (w+Aq - igBoY^-^-!^ - (w_Ao - 

r 

„ {w+Bo + z 9 Ho)e -**“'-/2 - {w_Bo A igHo)e -*‘“’+/2 

B = - , 

r 

where r = sjup' — and w± = {w ± ■r)/2. For 
\q/w\ <C 1 one has A ~ Aq, B ~ Bq, which is equiv¬ 
alent to dA/dt ~ 0, dB/dt ~ 0 in equation (fT2)l . pro¬ 
vided one neglects the “frequency shifts” of the order 
0(q^/w). This observation constitutes the method of 
RWA [sl, 113, which allows one to neglect the rapidly 
oscillating terms in differential equations for probability 
amplitudes. Hence, for /3n ^ 70 one can neglect the sec¬ 
ond and the third lines in equations (O - (ITTl) and the 
fourth line in equation (ED. Moreover, for nondiagonal 
elements one can also neglect the fourth line in equa¬ 
tions ® - (nni- Similar reasoning holds for the atomic 
and cavity dampings. 

Therefore, by choosing resonant modulation frequen¬ 
cies and performing RWA one obtains time-independent 
coupled differential equations for the matrix elements of 
p that can be solved analytically or numerically. The gen¬ 
eral solution can be quite cumbersome, so in this work 
we shall pursue analytically only the asymptotic behav¬ 
ior, obtained by setting the left-hand side of equations 
® - (HID to zero. 
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III. TWO-EXCITATIONS BEHAVIOR 


A. Photon generation from vacuum in resonant 
regime 


First we study the photon creation from the initial 
zero-excitation state (ZES) \g, 0). In the resonant regime, 
A_ = 0, the resonant modulation frequency was found 
two decades ago 0 and reads 

= 2wo + 7^goV^ , = ± . (14) 

The effective Hamiltonian is H = 6 \(pQ _){ip 2 ,n\ + h.c., 
where 

e = (15) 


^(r) Jr) (r) (r) 

= + -^+7^*72^. (16) 

2uJo 2 u!o go go 

Here denotes the dimensionless collective modula¬ 
tion depth [ 3 ^ . The analytical solution for this Hamilto¬ 
nian in absence of losses is straightforward [see equations 
(I551) - ((5711 below] and consists of periodic oscillations be¬ 
tween the states |</5o,-) ^ \V 2 ,n)- 

In the presence of dissipation described by SME, un¬ 
der the pure atomic damping we obtain the nonzero 
asymptotic probabilities of dressed states {pj^ = 
{‘fN,T\p\‘PN,T)) 


Po = 


i^r + (7/4)- 
3|0|^ + (7/4) 


2 ^ Pf = 


i^r 


3 | 0 |^ + ( 7 / 4 )" 


(17) 


and pf = \p^- For the pure atomic dephasing we obtain 
Po — P 2 — \ the pure cavity relaxation 


Po 


Pf 


5 |6»| V (3«;/4)^ ’ ^ 5|6>|V(3k/4)^ 

i(\/2 -f ifpl^ , = i(\/2 - ifpf . 


Hence dissipation affects dramatically the unitary dy¬ 
namics, although excitations can still be created when¬ 
ever K, 7 | 0 |. 

For our tentative parameters we obtain \ 6 \ ^ 10“^po; 
so the periodic generation of the state |<P 2 , 7 ?,) could 
be observed without difficulties in current experiments. 
This is illustrated in figure [T^ where we use the pa¬ 
rameters uJol^TT = 8 GHz, Po = 5 X 10“^wo, A_ = 0, 
en = 5 X 10“^Hoj = 2a;o -I- and consider mod¬ 

erate dissipation rates k = 7 = 70 = 2x lO'^^po- We 
plot the average photon number (h), the Mandel factor 
Q = [((Ah)^) —(h)]/(fi) (that quantifies the spread of the 
photon number distribution) and the atomic excitation 
probability Pe- We see that both models of dissipation 
predict identical results and several oscillations of (h), Q 
and Pe with high visibility could be observed within the 
time interval of 500 ns. 



FIG. 1: Photon generation from ZES in the resonant regime, 
(a) One-tone modulation fsection IlII Al) . (b) Two-tone mod¬ 
ulation (section II V Al) . Black (red) lines stand for SME 
(DPME). 


B. Anti-Jaynes-Cummings (AJC) regime 

Now we consider the photon generation from ZES in 
the dispersive regime, |A_|/2 ^ poV^max, where rimax 
is the maximum number of excitations. The resonant 
modulation frequency is 

= A+ - 2 ((5_ - S+) (19) 

where 6 - (5+) is the standard dispersive (Bloch-Siegert) 
shift given by (5± = go/A± and A+ = (^ -|- Oq) (we ne¬ 
glected higher order corrections to [^)- The unitary 
dynamics is described by the Anti-Jaynes-Cummings 
(AJC) Hamiltonian which in the dressed pic¬ 

ture reads H = 6 \ipo,v){'p 2 ,-v\ + h.c.^ 

e = i^goPT^-^^ ( 20 ) 


y(j) 


Ccj 





90 



( 21 ) 


Here V = A„/ |A_| = ± is the “detuning symbol” and 
we denote |(po,X)) = Iv^o.-)- 

For SME we find that under pure atomic dephasing the 
asymptotic solution is p^ = pf = with all other prob¬ 
abilities equal to zero. In the presence of pure atomic and 
cavity dampings we get the asymptotic nonzero proba¬ 
bilities (to the second order in go/A-) 


P 2 


-V 


^(l-#)+7# 

«|i+7(l-3||-) 


D 

Pi 


( 22 ) 






























5 


-V 

Pi 


^(l +#) 



(23) 


Po = 


' K + 7 + 2^ {k- l)' 

W\ 


P2 


-V 


(24) 


where Pq can be easily found from the normalization con¬ 
dition. Once again the dissipation changes drastically the 
dynamics, but for /c ,7 <C \9\ excitations are still gener¬ 
ated. In particular, for pure atomic (cavity) damping 
and 7 <C |P| (k <C |P|) we obtain asymptotically p^ ~ \ 
{Pi^ ~ 1 ), so only one photon (one atomic excitation) is 
generated. 

For our tentative parameters we obtain |0| ^ of 

the same order of magnitude as in section IIII A 1 so this 
effect could also be implemented in circuit QED with 
current technology. This is illustrated in figure [5^ for 
parameters a;o/27r = 8 GHz, go = 5 x = 8 g 0 j 

en = 5 X 10“^flo, = A 2 ,- — Aq -I- 26+ x 0.954 and 

K = 7 = 70 = go X 10“"^. Once again the results are 
almost indistinguishable for the two dissipation models 
and the periodic generation of one cavity and one atomic 
excitations is easily observable on timescales of the order 
of 500 ns. We notice that two excitations can also we 
generated from nonvacuum initial states for = A+ — 
2k(5- —15+), where k is an integer and the corresponding 
generation rate in equation (1^ becomes 9\/k [s^. 



FIG. 2: Photon generation from ZES in the dispersive regime, 
(a) AJC behavior ('section llllBIl . (b) DCE behavior (section 
0. Black (red) lines stand for SME (DPME). 


C. Anti-DCE regime 


The Anti-DCE behavior was predicted recently in [s^. 
It occurs in the dispersive regime and consists of coherent 
annihilation of two system excitations via the coupling 
|<Pfc,x>) t-)- \{pk- 2 ,-v) for a given value of fc > 3, roughly 
equivalent to the transition \g,k) -o- \e,k — 3). Eor the 
initial state |g)(g| ® Pfleid three photons are subtracted 
from the field due to the time-modulation of the system 
parameters, so the denomination Anti-DCE seems ap¬ 
propriate. Anti-DCE is implemented for the modulation 
frequency 

= 3wo - Oo -I- 2 (^_ - 5+) (fc - 1) , (25) 

where fc > 3 and the dressed-picture Hamiltonian is H = 
9*\'Pk,-D){Pk-2--D\ + h.C; 

e = ^Jk {k - 1) (fc - 2 )T^'^) ( 26 ) 

2cuq/\— 


Y(-4) 


£ 


(A) 


2 uj(} -f A_ 


,, , A JA) JA) 

^0 + _£g 

2a;o -l- A_ Dq go 


(27) 


This phenomenon occurs only for the modulation of w, 
H or g, so it cannot be implemented by the parametric 
down-conversion process in which only 0. For our 

tentative parameters, |A_| ^ lOgo and fc ^ 5, we get 
|0| ~ 10“^go, so Anti-DCE is a rather weak effect that 
requires fine-tuning of the modulation frequency and a 
prolonged maintenance of perturbation. Still, the tran¬ 
sition rate |61| is slightly greater than the state-of-the-art 
dissipation rates, so this effect lies on the threshold of 
implementability with current technology. 

For SME, under pure atomic dephasing we obtain 
asymptotically p^ = p ^_2 and = pf for 1 > 0 , so 
for large times the dephasing changes completely the ex¬ 
pected behavior. Anti-DCE behavior is also strongly af¬ 
fected by the cavity relaxation, when one gets approxi¬ 
mately 

(A + l)(^l-||-^g^+i (28) 

+ ^PNli-N(l-^y- 

= + (29) 

hold for the pure atomic damping: 
(l-|i(2A^+l))p)^?i (30) 

+ (-^ + P^+i ~ '^'^Pn 


dpN 

d{Kt) 


dpN,N 

d{Kt) 


Similar formulae 

dp% ^ 
d{lt) 
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d{lt) 



-V 

N-\-l 



— V 

Pn ■ 


(31) 


In both cases the system state goes asymptotically to 
ZES. This did not happen in the phenomena analyzed in 
sections IlII Al and IIII Bl because there the ZES was cou¬ 
pled to 2 -excitations states by the modulation, resulting 
in nonvacuum equilibrium state. 



0.6 
Q 0.4 
0.2 
0.0 

pO-3 

eo .2 

0.1 

0.0 



t(|JS) 


t(ps) 


FIG. 3: Coherent photon annihilation in the dispersive regime 
for initial coherent state \g) 0 |a). (a) Anti-DCE behav¬ 

ior fsection IlII C|l . (b) Enhanced Anti-DCE behavior (section 
HVBl) . Black (red) lines stand for SME (DPME). 


Figure |3^ assesses the feasibility of experimental obser¬ 
vation of Anti-DCE for the initial coherent state jg)® ja). 
We use the values a = 2, ujqI2tt = 8 GHz, go = 
5 X lO^^wo, A_ = 8 go, en = 5 x 10“^Oo, k = 10“®go, 
7 = 7 ,^ = 5 X 10“®go and 774 "^^ = A 4 _+ — A 2 ,- — 2(5+ x 2.791. 
Both master equations predict practically identical re¬ 
sults and for initial times the system exhibits the ex¬ 
pected behavior, characterized by the periodic reduction 
of the average photon number accompanied by the partial 
excitation of the atom and modification of the field statis¬ 
tics from Poissonian to super-Poissonian. The curves cor¬ 
responding to (n), Q and Pe exhibit broad widths due to 
fast oscillations associated with off-resonant exchange of 
excitations between the atom and the field that occurs on 
timescales (~ |A_|“^) much smaller that the character¬ 
istic time of Anti-DCE. For larger times the system tends 
to the ZES and the amplitude of oscillations, as well as 
the associated visibility, decrease. Yet the main features 
of Anti-DCE could be resolved in experiments with the 
current state-of-the-art parameters on timescales of a few 

/iS. 


IV. FOUR-EXCITATIONS BEHAVIOR 

A. Enhanced photon generation from vacunm in 
resonant regime 

Now we consider the two-tone multi-modulation ap¬ 
plied simultaneously to H and g, with one frequency given 
by equation (HI and the second 

= 2a;o -f goV2{V2n2 -'ll) , IR., 7^2 = ± . (32) 

For initial ZES the effective Hamiltonian reads H = 
0 \(po-){‘P2,Ti\ + d2\^2,Ti){‘f4,n2 \ + d.c., where 9 is given 
by (fTSl) and 

d2 = 7go7e2^T(’'2) (33) 

(r2) (r2) (r2) (r2) 

X(r2) = ^ ^ + 77e2^(2 + 'R'R 2 V 2 ) . 

2uJq 2lOq Qq Qq 

(34) 

It is worth writing first the solution for the unitary 
evolution. We expand the interaction-picture wavefunc- 
tion as \'ip) = A\g)Q_) + B\ip 2 ,n) + C|7^4,TC2) + E(- ’ 0. 
where X)(''') denotes the other terms not coupled by the 
Hamiltonian and we assume real probability amplitudes. 
For the initial condition Hq = 0 we find 

41 = (dlo + cos Rt + ^ [29*e*Aa (35) 

-2d*(r Ao + 02 CO) cos^ ^Rt - (|0|2 - |d2|2)Co 


B = -i 


{9*Ao + 92 C 0 ) sin Rt 
R 


(36) 


C = 


29 ^( 9 *Ao +6200) cos^ 
m^-\d2\^)Co-‘29*9*Ao 


(37) 


where R = \/\9\^~+~\9^. One can transfer the 
populations entirely from {A,C} to C at the small¬ 
est time Rtmin = TT by choosing 02 = x9*, x = 
(^Cq ± ^/Cq + Aq^ /Aq. In the present case x = ±1, so 
one can create the pure dressed state at the time 

tmin = 'K/{\9\y/2). 

For the pure atomic dephasing described by SME we 
get in the asymptotic limit = pf = 1/5, while 

for the pure atomic damping the nonzero probabilities 
are 


( 2 \ 9 \^ + (7/4)^ _ , 101 ^-(7/4)^ \ 

1 , | 0 |^ 3|0|^+7V4; 


pf (38) 
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Pi — 


m 


3 | 6 >r+72/4 




(39) 


+ _ 1 5|6>|V72/4 ^^ 
2 3 | 0 | 2 + 72 / 4 ^^ 


(40) 


r 3 ^2 2 


(41) 


where can be found from the normalization condi¬ 
tion. Similar expressions are obtained for the pure cavity 
damping, so we conclude that for k ,7 <C | 0 | four excita¬ 
tions can be generated from vacuum. For our tentative 
parameters we get | 0 |,| 6 > 2 | 10 “^ 5 o, so this effect can 

also be implemented with current technology. 

We illustrate the typical outcome in figure [TJ), where 
we consider the modulation of fl with frequency 77 ^’’^ = 
2 u!o + go \/2 and modulation of g with frequency = 
2a;o + go{2 — \/2). Other parameters are wq/Stt = 8 GHz, 
go = 5 X lO^^wo, A_ = 0, eo = 5 X 10“^Ooj £g = 1-97 x 
= 1 and K = 'y = 'y^ = 2 x 10 “"'’go- 
We see that even for moderate values of dissipation rates 
the enhanced photon generation due to multi-modulation 
can still be observed. Moreover, the predictions of SME 
and DPME are almost indistinguishable. 


to at specific times without affecting the other 

states, so four photons are indeed annihilated. If the ini¬ 
tial statistics is not known, the transfer is not complete 
because the required value of O 2 cannot be determined 
beforehand. 

For SME we reach the same conclusions about the 
asymptotic behavior as in section fill Cl supplemented by 
the additional condition p® = p ®_2 = P ^-4 under pure 
atomic dephasing. As the associated transition rates are 
similar to the one estimated in section ITlI Cl the enhanced 
Anti-DCE also lies at the threshold of experimental feasi¬ 
bility, provided two simultaneous fine-tuned modulations 
can be sustained for a time interval of a few ps. This 
is illustrated in figure [3b for the initial coherent state 
|g) C) I a) and the simultaneous modulation of H and g. 
We used the parameters of figureEb, Sg = 9.91 x 10“'^go, 
^(4i2) _ ^ — Aq -I- 25+ X 0.954 and nonzero modulation 

weights = 1. We see that the reduction of 

(h) is stronger than in the standard Anti-DCE case, oc¬ 
curring at larger timescales; the Q—factor increases while 
the atomic excitation probability undergoes faster oscil¬ 
lations and is slightly smaller than in figure |3b because 
the transfer of population from |g, 4) to |e, 1) is only par¬ 
tial. Moreover, SME and DPME predict almost identical 
behaviors, with minor quantitative differences appearing 
only for larger times t ^ 5 ps. 


B. Photon annihilation in enhanced Anti-DCE 
regime 


V. DYNAMICAL CASIMIR EFFECT (DCE) 
WITH SINGLE QUBIT 


One can enhance the Anti-DCE behavior by combining 
the modulation frequency (j25p with the second frequency 

=A+-2(5_-5+)(A:-3) , (42) 


In the dispersive regime one can create many pairs of 
excitations from the initial ZES for the modulation fre¬ 
quency 

=2(a;o + <5--<5+). (45) 


where fc > 4 is an integer. The effective Hamiltonian be¬ 
comes H = 9 *\ipk.v){g}k- 2 -v \ + 9 l\ipk- 2 -v){Vk-i,T' \ + 
h.c., where 9 is given by (IMl) and 

92 = i^goVVk^T^^^^ (43) 


Y(-42) _ 


(A2) (A2) {A2) 

_|_ ^9_ 

A+ A+ go 



(44) 


In this regime one combines simultaneously the cou¬ 
pling \(pk,Ti) —> \^k- 2 ,-Ti) via Anti-DCE with the cou¬ 
pling \(pk- 2 ,-v) —^ l^k-i,!^) via the AJC behavior, ob¬ 
taining the annihilation of four excitations. In the disper¬ 
sive regime this corresponds to subtraction of four pho¬ 
tons from the state |g, k){g, k\. The effective Hamiltonian 
is similar to the one considered in section ITV A[ so the for¬ 
mulae (I 35 ]) - (l37l) can be applied after appropriate substi¬ 
tutions. If the initial cavity state is known and the pop¬ 
ulation of the excited atomic state is negligible, one can 
transfer completely the populations of {|<Pfc,r>), |<Pfc- 4 ,x>)} 


The dynamics is described by the effective Hamiltonian 

00 

771=0 

(46) 

^ (47) 


Jd) (d) (d) . (d) 

ojQ Do go Do 5- 


(48) 


so only the states \ 1 p 2 m ,T>) ^ |g, 2 m) are populated during 
the dissipationless evolution. The unitary dynamics of 
the Hamiltonian was studied in , where it was shown 
that the photon generation from vacuum suffers satu¬ 
ration due to effective Kerr nonlinearity (of magnitude 
gl/Al), and ( fi) exhibits a collapse-revival behavior as 
function of time. Simple analytical expressions cannot 
be obtained for this case because many dressed states 
















are coupled simultaneously and the argument of the ex¬ 
ponentials depends on the index m in nontrivial manner. 

For SME the pure atomic damping is equivalent to 
the pure cavity damping with effective relaxation rate 
Kef = 7 ( 50 /^-)^ (for the initial ZES and under RWA). 
The asymptotic solution for a similar problem (paramet¬ 
ric amplification in the presence of Kerr nonlinearity and 
cavity relaxation) was obtained two decades ago using 
the method of potential solutions for the corresponding 
Eokker-Planck equation [d^. However, the asymp¬ 
totic solution is of little use in our case because the pho¬ 
ton statistics during the transient time can be quite dif¬ 
ferent from the asymptotic one [s^ . so we rely entirely 
on numerical simulations. 

For our tentative parameters we estimate ^ 

so in principle the single-qubit DCE could be ob¬ 
served in the current state-of-the-art architectures. This 
is confirmed in figure [5 ]d where we show the expected be¬ 
havior for parameters a;o/27r = 8 GHz, go = 5 x IO^^wq, 
A_ = 8go, en = 5 X 10“^Ho, k = 7 = 70 = 5x 10“®go 
and = A 2 ,+ — Aq — 2(5+ x 1.02. Both dissipa¬ 
tion models predict identical results for initial times; for 
times t ^ 1 /iS some small quantitative differences appear 
without changing the overall behavior. The saturation 
in photon growth occurs due to the nonzero values of 
{^m,+2,V — Am. 15 — fjid)) 

oc m in (H51) for large values of 
m, and is not related to the presence dissipation. We see 
that photons can be generated and the collapse-revival 
of (h) can be observed on the timescale of a few /is. No¬ 
tice that the collapse-revival behavior is not associated 
to the absorption and reemission of photons by the atom 
[U , as the atomic excitation probability stays < 8% 
for all times. More photons can be generated if one adds 
a second modulation frequency or increases or A_ 

[ 3 ^ . but in the latter case the parameter \9m\ decreases 
so the process becomes slower and demands more precise 
tuning of the modulation frequency. 

VI. CONCLUSIONS 

We estimated the rates of photon creation and anni¬ 
hilation for some nonstationary phenomena induced by 


fast modulation of qubit parameters in dissipative circuit 
QED. Two different Markovian master equations (“stan¬ 
dard” and “dressed-picture” master equations) were used 
to account for the common sources of dissipation, and we 
verified that for the relevant regime of parameters the 
predicted behaviors are almost identical. 

Our results indicate that all the analyzed phenomena 
could be implemented experimentally with the present 
technology, provided one can implement modulations 
with relative depths |T| ^ 10“^ and frequencies rj ^ 
15 GHz. The most accessible phenomena (“fast phenom¬ 
ena”) are the two- and four-excitations generation from 
vacuum in the resonant regime (using single- and two- 
tone modulations, respectively) and the two-excitation 
generation via the AJC behavior in the dispersive regime. 
For these effects the photon generation rates are |0| ^ 
10“^goj so they could be implemented even in systems 
with moderate dissipation. On the other hand, genera¬ 
tion of several photons from vacuum via DCE and annihi¬ 
lation of two or four excitations via Anti-DCE in the dis¬ 
persive regime (“slow phenomena”) occur at lower rates 
|0| ^ so only state-of-the-art architectures would 

allow for their experimental verification. Additionally, 
the modulation frequency must be fine-tuned with ac¬ 
curacy of the order of |0| and maintained for the time 
interval ^ 500 ns for the fast phenomena and ~ 5 /its for 
the slow ones. Although all these requirements are chal¬ 
lenging, the verification of the nonstationary phenomena 
in circuit QED could lead to new manners of generating 
entangled states by means of the counter-rotating terms 
in the light-matter interaction Hamiltonian. 
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